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ABSTRACT 

Formulae are presented for the linear growth factor D/a and its logarithmic deriva- 
tive dlnD/dlna in expanding Friedmann- Robertson- Walker Universes with arbitrary 
matter and vacuum densities. The formulae permit rapid and stable numerical evalua- 
tion. A fortran program is available at http://casa.colorado.edu/~ajsh/growl/. 
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1 INTRODUCTION 

The linear growth factor g = D/a, where D is the amplitude of the growing mode and a is the cosmic scale factor, determines 
the normalization of the amplitude of fluctuations in Large Scale Structure (LSS) relative to those in the Cosmic Microwave 
Background (CMB) (Eisenstein, Hu & Tegmark 1999 Appendix B.2.2). Its logarithmic derivative, the dimensionless linear 
growth rate / = dlnD/dlna, determines the amplitude of peculiar velocity flows and redshift distortions (Peebles 1980 §14; 
Willick 2000; Hamilton 1998). As such, the growth factor g and growth rate / are of basic importance in connecting theory 
and observations of LSS and the CMB. 

In a Friedmann-Robertson- Walker (FRW) Universe containing only matter and vacuum energy, with densities f2 m and 
£Ia relative to the critical density, the linear growth factor is given by* (Heath 1977; Peebles 1980 §10) 

where a is the cosmic scale factor normalized to unity at the epoch of interest, H (a) is the Hubble parameter normalized to 
unity at a — 1, 

H(a) = (fi m a~ 3 + tt k a- 2 + f2 A ) 1/2 , (2) 
and the curvature density Q k is defined to be the density deficit 

fik = 1 — f2 m — Qa , (3) 

which is respectively positive, zero, and negative in open, flat, and closed Universes. The normalization factor of 5fi m /2 in 
equation (1) ensures that g — *• 1 as a — *• 0. It follows from equation (1) that the dimensionless linear growth rate / is related 
to the growth factor g by 

The fact that the integrand on the right hand side of equation (1) is a rational function of the square root of a cubic implies 
that the integral can be written analytically in terms of elliptic functions. However, the analytic expressions are complicated, 

* Regarded as a function of cosmic scale factor a, the growth factor evolves as 
5f2 m H(a) f" da' 5!! m (a)« 2 /f(«) 3 f a da' 

g{a) - 



2 a J a' 3 H(a') 3 2 J a' 3 H{a') 3 

where Q m (a) = S7 m a _3 /_ff(a) 2 . 
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and have yet to be implemented in any published code. Explicit analytic expressions in the special case of a flat Universe, 
fik = 0, have been given in terms of the incomplete Beta function by Bildhauer, Buchert & Kasai (1992) and in terms of the 
hypergeometric function by Matsubara (1995). Analytic expressions for the luminosity distance in the general non-flat case 
are given in terms of elliptic functions by Kantowski, Kao & Thomas (2000). 

Lahav et al. (1991) give a simple and widely used approximation to the growth rate^: 

/(n m , n A ) « Q.T + (1 + On/2) Oa/70 . (5) 

From this, together with relation (4), follows the approximate expression for the growth factor quoted by Carroll, Press and 
Turner (1992): 

«?(On, n A ) « — T . (6) 



Q.T - Q A + (1 + On/2)(l + Oa/70) 

Given the increasing precision of measurements of fluctuations in the CMB (de Bernardis et al. 2000; Hanany et al. 2000) 
and LSS (Gunn & Weinberg 1995; York et al. 2000; Colless 2000) and the growing evidence favouring a cosmological constant 
(Gunn & Tinsley 1975; Perlmutter et al. 1999; Riess et al. 1998; Kirshner 1999), it seems timely to present exact expressions, 
suitable for numerical evaluation, for the growth factor g (hence /, through formula [4]) valid for arbitrary values of the 
cosmological densities Q m in matter and Qa in vacuum. 

The procedure presented in this paper is to expand the integral of equation (1) as a convergent series of incomplete Beta 
functions, conventionally defined by 



B(x,a,b) = f f-^l-t) 6 - 1 ^ . 
Jo 



(7) 



In effect, the method can be regarded as generalizing Bildhauer et al.'s (1992) formula. What makes the scheme attractive is 
that the Beta functions in successive terms of the series can be evaluated recursively from each other through the recursion 
relations 

aB(x,a,b) = x a (l - x) b + {a + b) B(x, a + 1, b) (8) 

bB{x,a,b) = -x a (l-x) b + (a + b)B(x,a,b + l) (9) 

aB(x,a,b+l) = x a (l- x) b + bB(x,a+l,b) (10) 

the last of which follows from the other two. In practice, each evaluation of the growth factor g involves from one to seven 
calls to an incomplete Beta function, followed by elementary recursive operations. 

The resulting numerical algorithm is fast (provided that 57 m and \£Ia\ are not huge - see §3) and stable, and has been 
implemented in a fortran package growA available at http : //casa. Colorado . edu/~ajsh/growl/. The growA package includes 
an updated version of an incomplete Beta function originally written by the author a decade ago. 



2 FORMULAE 

Figure 1 shows contour plots of the growth factor g(ftm, Qa) and growth rate /(fim, Qa) computed from the formulae presented 
below, as implemented in the code growA. The results have been checked against numerical integrations with the Mathematica 
program. 

Perhaps the most striking aspect of these plots, emphasized by Lahav et al. (1991), is that the growth rate / is sensitive 
mainly to Q m , with only a weak dependence on Qa- 

Figure 2 shows the ratio of the Lahav et al. (1991) growth rate /, equation (5), to the true growth rate. The Figure 
illustrates that the Lahav et al. approximation works well except at small Q m . The approximation (with the 4/7 exponent 
advocated by Lightman & Schechter 1990) works particularly well for currently favoured cosmologies, being accurate to better 
than 1% for flat Universes with f2 m = 0.2-3.9. 



2.1 Limiting cases 

The requirement that the Hubble parameter H(a), equation (2), be the square root of a positive quantity for all a from 
(Big Bang) to 1 (now) imposes two requirements. The first is that the matter density be positive 

t Lahav et al. adopt Q 0,6 rather than f2 4 / 7 , following Peebles (1980). The 0.6 exponent is more accurate for low f2 m , but 4/7 (Lightman 
& Schechter 1990) works better elsewhere, and in particular is more accurate for currently favoured cosmologies, flat Universes with 
ilm > 0.2. The 4/7 exponent is therefore currently favoured. 
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Figure 1. Contour plots of (left) the growth factor g(Q m , Ha), and (right) the growth rate /(C m , Ha), in expanding Universes containing 
matter and cosmological constant with densities f! m and Oa relative to the critical density. Universes below the 45° dashed line are 
geometrically open, while those above are closed. Universes to the upper left of the long dashed line are decelerating, while those to 
lower right are accelerating. Universes to the left of the almost vertical line near SIa ~ will eventually turn around, while Universes to 
the right will expand forever. The region at the bottom right corner is physically inaccessible to Universes that expand from zero and 
that contain only matter and cosmological constant. The boundary between turnaround and eternal expansion, and the boundary of the 
inaccessible region, together form, in the closed case, the approaching and receding parts of the loiter line. Starting at Sl m = 1, Ha = 0, 
a Universe can evolve upward and rightward along the approaching part of the loiter line to Einstein's loiter point at f2 m = 2 Ha — ► oo. 
After an indefinite period of hanging around, the Universe can then either recollapse along the same loiter line, or else continue into 
renewed expansion downward and leftward along the receding part of the loiter line, the boundary of the inaccessible region. The growth 
factor g tends to infinity at the boundary of the inaccessible region, but only contours up to g = 5 are marked. The contour plot of 
g(Qm,Qj^) is similar to Figure 7 of Carroll et al. (1992). 




Figure 2. Contour plot of the ratio of the Lahav et al. (1991) growth rate /, equation (5), to the true growth rate. The plot illustrates 
that the Lahav et al. approximation works well except for small f2 m . In particular, the approximation is accurate to better than 1% 
for currently favoured cosmologies, flat Universes (45° dashed line) with Q m > 0.2. For small S7 m (and any Ha) a somewhat better 
approximation results if the Lightman & Schechter (1990) exponent 4/7 in equation (5) is replaced with 0.6, as in Lahav et al.'s original 
paper. The white region in the diagram at Q m <J 0.1 and Ha < 1 is where the error in the Lahav et al. approximation (with the 4/7 
exponent) exceeds 10%; contours in this region are omitted to avoid confusion. 
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fi m > o , 

and the second can be interpreted as a condition that the Universe not be too closed 



3Q n 



1/3 



(11) 



(12) 



If either of the two conditions (11) or (12) were violated, then the Hubble parameter H(a) would go to zero at some finite 
cosmic scale factor a, indicating that the Universe did not expand from a — 0, but turned around from a collapsing to an 
expanding phase at some finite a. 

The limiting values of the growth factor g and growth rate / as the matter density goes zero, 57 m — > 0, are 



g (Q m = o,q a ) = , /(n m = o,n A ) = o 



(13) 



provided also that Qa < 1, in accordance with the condition f2 k > from equation (12). Physically, structure cannot grow in 
a Universe containing only vacuum. 

The second condition, equation (12), is saturated when Q^ — ► — (27Q 2 i1 Qa/4:) 1 ^ 3 with Q^ < — 3 m /2. This marks the 
boundary of the inaccessible region to the bottom right of Figure 1. The limiting values of the growth factor g and growth 
rate / in this case are 



g(Q m , Qa) — » oo , /(fi m , Qa) 



(14) 



Physically, the growth factor tends to infinity because such Universes are in renewed expansion after having spent an indefinite 
period of time at Einstein's unstable loitering point, at O m = 2Qa — > oo. 



2.2 Case |S7 k /(l - fi k )| < 1 

For small curvature density fik, expand the integral in equation (1) as a power series in fik: 
n(Q SOn f da (-)"(3/2) w / » k a- 2 V 

g(^,M A )- 2 ^ a 3 (fima -3 + nA) 3/ 2 Z. „. ^fi m a-3 + n A J (15) 

where (x)„ = x(a; + l)...(a; + ri — 1) is a Pochhammer symbol. 

Each term of the sum on the right hand side of equation (15) integrates to an incomplete Beta function. The cases of 
positive and negative cosmological constant must be distinguished. For positive cosmological constant, Qa > 0, 



while for negative cosmological constant, f2 A < 0, 



Q m + fi A ' 6 3 ' 3 3 



5fi^V^ (-)"( 3 /2)n Z' ^ V D / 5 , n 2 , 2n 

'6 t 
i 

CO! 

5nL /3 y> (-)"(3/2)„ / O k \" 5 n _1 



+ ? ,tt + ^ , (16) 



In each of formulae (16) and (17), incomplete Beta functions must be evaluated for three terms, and then the remaining Beta 
functions can be evaluated recursively from these three, through the recursion relations (8)-(10). In equation (16) (Qa > 0), 
the recursion is stable from n — upward, or from large n downward, according to whether the argument Q A /(f2 m + Qa) is 
greater than or less than"!" 1/3 j n equation (17) (Qa < 0), the recursion is stable from n — upward or large n downward as 
the argument |f2 A |/fi m is greater or less than 1 + (3/2) [(>/2 - 1) 1/3 - (\/2 - 1)~ 1/3 ] = 0.106. 

How fast do the expansions (16) and (17) converge? Convergence is determined essentially by the convergence of the parent 
expression (15) at the place where the expansion variable Q^a^ 2 / (f2 m a~ 3 + Qa) attains its largest absolute magnitude over the 
integration range a G [0, 1]. For Qa > 0, the expansion variable attains its largest magnitude at a = min{l, [f2 m /(2 O/v)] 1 / 3 }, 
while for Qa < 0, the expansion variable is always largest at a — 1. Physically, the extremum at f2 m a -3 = 2fi A occurs where 
the Universe transitions from decelerating to accelerating. It follows that if Q m > 2Qa (decelerating), then successive terms 
of the expansions (16) and (17) decrease by w Q ii /(Q rn + Qa), whereas if Q m < 2 Qa (accelerating), then successive terms 

t A somewhat lengthy calculation, confirmed by numerical experiment, shows that the asymptotic (i.e. after many iterations) point of 
neutral stability of the recurrence B(x, a, b) — > B(x, a + m,a + n), where m and n are positive or negative integers, occurs at that unique 
point x € [0, 1] satisfying 

m m n n 

x m (l-x) n = ± - " . 

(m + n) m + n 
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decrease by « 2 2/3 fi k /(3 fi 2 / 3 !^ 73 ). Thus n terms of the expansions (16) and (17) will yield a precision of w [nk/(fi m + fiA)] r ' 
if fi m > 2fi A , or a precision of « [2 2/3 fi k /(3 fi^fi^ 3 )]™ if !i m < 2fi A . 

2.3 Case |fi A /(l - fi A )| < 1 

For small cosmological constant fiA, expand the integral in equation (1) as a power series in Oa: 
5fi m f 1 da <A H n (3/2)„ / fi A 



(O ^- 5 " m / ™ \^ (~) (3/2)n ^ »A V 

51 m ' Aj ~ 2 A, ^(Q^ + ^^pL n! U m a-3 + n k a- 2 J • 

J0 V n=0 

Again, each term of the sum on the right hand side of equation (18) integrates to an incomplete Beta function. The cases 
of positive and negative curvature density fi k must be distinguished. For an open Universe, fi k > 0, 

n(Q o _ 5»jL y> (~)"(3/2)„ /n^i A V / _fik 5 \ 

while for a closed Universe, fi k < 0, 

5fi?n f H"(3/2)„/^aV Jiai 5 |9 1 



n=0 



9(fim ' fiA) = 2N^L — ^ — [jkF B \^:'2 +3n '-2- n ■ (20) 



Here the Beta functions in successive terms can be computed recursively from a single Beta function. In equation (19) (fi k > 0), 
the recursion is stable from n = upward or large n downward as the argument fi k /(fi m + fi k ) is greater than or less than 
3/4. In equation (20) (fi k < 0), the recursion is stable from n — upward or large n downward as the argument |fi k |/fi m is 
greater than or less than (3/2) [(a/2 + 1) 1/3 - (\/2 + 1)~ 1/3 ] = 0.894. 

The convergence of the expansions (19) and (20) is determined by the convergence of the parent expansion (18) at the 
place where the expansion variable fi A /(fi m a + fi k a~ ) attains its largest magnitude over the integration range a £ [0, 1], 
which occurs at a = 1 for both positive and negative fi k . It follows that successive terms of the expansions (19) and (20) 
decrease by w fi A /(fi m + fik), and n terms will yield a precision of « [fi A /(fi m + fi k )] n . 



2.4 Case |fi m /(l - fi m )| < 1 



For small matter density fi m , a power series expansion of the integral in equation (1) is again possible, but the integral must 
be split into two parts to ensure convergence of the integrand over the full range a £ [0, 1] of the integration variable. Define 
the cosmic scale factor o eq at matter-curvature 'equality' to be where |fi m a<^ 3 /(fi k a^q 2 + Ha)| = |fik£i^iq 2 /(^ m0 eq 3 + fiA)|, and 



similarly that at matter-vacuum 'equality' to be where |fi m o eq /(fi k o eq + fi A )| = |fiA/(fi m a cq + fika eq 



These conditions 



d c q fi k tt c q 



if fi k > 0, and fi m a cq = fiA if fiA > 0. A good procedure is to use one of the methods of the 
§2.2 or §2.3, to integrate up to the cosmic scale factor a oq at matter-curvature or matter-vacuum 



reduce to fi 

previous two subsections 
'equality', whichever is later (larger a cq ) since the later epoch yields a more convergent fi n 
integral using the small fi m expansion: 



g>(fi m , fi A ) = G(a cq ) + 



5fi n 



da 



a 3 (fi k a- 2 + fi A ) 3 / 2 



E 



R n (3/2)„ 



fin 



fi k O -2 + fiA 



series, and then to complete the 



(21) 



The first term on the right hand side of equation (21) is the integral up to a cq , evaluated by one of the methods of the previous 
two subsections: 



G(a eq ) = 



5fi m 

2 Jo 



da 



a 3 (fi m a- 3 + fi k a- 2 + fi A ) 3 / 2 H(a eq ) 



g [fi m (Oeq), fiA(a eq )] 



(22) 



with fim(a) = fi m a~ 3 /-ff(a) 2 , fik(a) = fik«~ 2 / H (a) 2 , fi A (a) = fiA / H (a) 2 , and H(a) given by equation (2). 

As in the previous two subsections, each term of the sum in the second term on the right hand side of equation (21) 
integrates to an incomplete Beta function. The cases of positive and negative curvature and vacuum densities must be 
distinguished. For an open Universe with a positive cosmological constant, fi k > and fiA > 0, 



ff(fi m , fiA) = G(a cq ) + 



5 fin 



40, o 1 / 2 

*"k" A n=0 V "k 

while for an open Universe with a negative cosmological constant, fi k > and fiA < 0, 



E 



(-) n (3/2) n /fi m fi 



,1/2 



fi 



3/2 



B 



fi k a 



fi k a -2 + fiA 



,1 + ' 



3n 



g(fi m ,fiA) = G(a oq ) + 



5 fin 



(-)"(3/2) n / fi m |fi A 



4fi k |fi A !V2 

n— 



|l/2 



fi 



3/2 



B 



[fiAj 

fi k a~ 



. n 

2 2' 2 



(23) 



(24) 
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The expression for a closed Universe with positive cosmological constant, fi k < and Qa > 0, turns out never to be useful, 
but for reference it is 

s( n m , n A ) - G(a cq) + _^ ^ — j , i + _,-__ „j 

The fourth option, a closed Universe with negative cosmological constant, f2 k < and Qa < 0, does not yield a convergent 
expansion in f2 m . The small 57 m expressions (23)-(25) are more complicated than the small f2 k and small Q.a expressions 
obtained in the previous two subsections, since equations (23)-(25) each involve a sum of not one but three expansions, one 
to evaluate the first term on the right hand side, one to evaluate the second term at a — a eq , and a third to evaluate the 
second term at a = 1. 

It will now be argued that a small f2 m expansion is advantageous only in the case where matter-vacuum equality occurs 
after matter-curvature equality. In particular, this has the consequence that the expansion (25) is never useful. Suppose that 
matter-curvature equality, |fi m a7 q 3 /(nkaiq 2 + ^a) = |f2 k a~ q 2 / (fi m a~ q 3 + Oa)|, occurs after matter- vacuum equality. Then the 
first term in equation (21), G(a eq ), would be evaluated by the small f2 k method, equation (16) or (17). However, if matter- 
curvature equality occurs after matter-vacuum equality, then it is also true that matter-curvature equality occurs after (larger 
cosmic scale factor a) the transition, at fi m a~ 3 = 2f2 A , from a decelerating to an accelerating Universe. As discussed in the 
last paragraph of §2.2, the convergence of the small fi k expansions (16) and (17) are then determined by the convergence 
of the parent equation (15) at the deceleration-acceleration transition, not by its convergence at a later epoch. Thus, if 
matter-curvature equality occurs after matter-vacuum equality, then there is no point in splitting the integral at a eq as in 
equation (21); one might as well use the small fl k expressions (16) or (17) all the way to a — 1, since they converge just as 
fast at a = 1 as at a = a eq . In fact the small fi k expressions (16) and (17) are computationally faster than the small fi m 
expansions (23)-(25), because the former involve a single sum whereas the latter involve three. 

The conclusion from the previous paragraph is that a small f2 m expansion is advantageous only in the case where matter- 
vacuum equality occurs after matter-curvature equality. Examination of evolution in the f2 m -fiA plane reveals that matter- 
vacuum equality happens after matter-curvature equality only if fi k > (conversely, matter-curvature equality happens after 
matter- vacuum equality only if S7a > 0). Thus, in the cases where the small O m expansion is useful, the relevant expressions 
to use are (23) or (24) with the first term, G(a eq ), being evaluated by the small 51a expansion (19) with fi k > 0. 

The Beta functions in successive terms of the expansions in the second term on the right hand sides of equations (23)-(25) 
can be computed recursively from two Beta functions. In equation (23) ($l k > and Oa > 0), the recursion is stable from 
n = upward or large n downward as the argument f2 k a~ 2 /(fi k a~ 2 + Oa) (with a = a oq or 1) is greater than or less than 
(3/2)[(-\/2 + 1) 1/3 - (a/2+ 1)~ 1/3 ] = 0.894, the same as for equation (20). In equation (24) (f2 k > and Q. A < 0), the recursion 
is stable from n = upward in all cases. In equation (25) (fi k < and Q,a > 0), the recursion is stable from n = upward or 
large n downward as the argument |f2 k |a~ 2 /r2A is greater than or less than 3/4, the same as for equation (19). 

The convergence of the expansions (23)-(25) is determined by the convergence of the parent expansion (21) at the place 
where the expansion variable O m a -3 /(fi k a -2 + Qa) attains its largest magnitude over the integration range a £ [o eq , 1], which 
occurs at a — a oq in all cases. It follows that successive terms in the expansions (23)-(25) decrease by w O m a ( I q 3 /(f2 k a^ q 2 + f2A), 
and n terms will yield a precision of w [f2 m aj q 3 /(n k aj q 2 + Aa)]" '■ 

2.5 Which formula to use? 

The various formulae (16), (17), (19), (20), and (23)-(25) converge over overlapping ranges of fi m and JIa- A sensible strategy 
would be to choose the expression that converges most rapidly. 

If fi k < (closed Universe), then use the small f2 k (§2.2) or small J7a (§2.3) methods as |O k /(l — Q k )| or |f2 A /(l — Qa)| 
is smallest. 

If fi k > (open Universe), then use the small f2 k (§2.2), small Qa (§2.3), or small f2 m (§2.4) methods as |fi k /(l — O k )|, 
|fi A /(l — or |n m /(l — n m )| is smallest. If |fl m /(l — fl m )| is smallest, then determine which occurs later (larger a eq ), matter- 
curvature equality |n m a7 q 3 /(r2 k a^ q 2 + Q A )| = |r2 k a~ q 2 /(f2 m a~ q 3 + S1a)|, or matter-vacuum equality |n m aj q 3 /(n k a^q 2 + Oa)| = 
\ClA/(0, m a~q + J7 k a^q 2 )|. If the former, then revert to using the small f2 k method, §2.2. If the latter, then use one of the small 
Qm expressions (23) or (24) with the first term, G(a eq ), equation (22), evaluated by the small Qa expression (19). 



3 COLLAPSING UNIVERSE 

The procedure proposed in this paper fails for Universes that are collapsing rather than expanding. As a Universe approaches 
turnaround, the Hubble parameter, and consequently the critical density, approaches zero, causing one or more of Q m , Q k , 
and r^A to tend to ±oo. In such cases the expansions presented herein converge ever more slowly. 



(25) 
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It is not clear how to adapt the method to deal with Universes near turnaround and thereafter, or indeed whether this 
is possible. 

4 SUMMARY 

Formulae suitable for numerical evaluation have been presented for the linear growth factor <j(f2 m ,fiA) = D/a and its 
logarithmic derivative /(f2 m ,f2A) = dlnD/dlna in expanding FRW Universes with arbitrary matter and vacuum densities 
f2 m and Q.a- A fortran package growA implementing these formulae is available at http: //casa. Colorado . edu/~ajsh/growl/. 
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